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Abstract 

We consider a class of dynamical systems on a compact Lie group G with 
a left-invariant metric and right-invariant nonholonomic constraints (so called 
LR systems) and show that, under a generic condition on the constraints, such 
systems can be regarded as generalized Chaplygin systems on the principle 
bundle G ^ Q — G/H, H being a Lie subgroup. In contrast to generic 
Chaplygin systems, the reductions of our LR systems onto the homogeneous 
space Q always possess an invariant measure. 

We study the case G = SO{n), when LR systems are multidimensional gen- 
eralizations of the Veselova problem of a nonholonomic rigid body motion which 
admit a reduction to systems with an invariant measure on the (co)tangent bun- 
dle of Stiefel varieties V{k, n) as the corresponding homogeneous spaces. 

For k = 1 and a special choice of the left-invariant metric on SO{n), we 
prove that after a time substitution, the reduced system becomes an integrable 
Hamiltonian system describing a geodesic flow on the unit sphere S"~^. This 
provides a first example of a nonholonomic system with more than two degrees 
of freedom for which the celebrated Chaplygin reducibility theorem is applicable 
for any dimension. In this case we also explicitly reconstruct the motion on the 
group SO{n). 

*AMS Subject Classification 37J60, 37J35, 70H06, 70H45 
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1 Introduction 

In classical nonholonomic mechanics a special attention is given to Chaplygin 
systems whose Lagrangian and constraints admit symmetries such that after an 
appropriate reduction the equations of motion take the form of unconstrained 
Lagrangian systems with some extra (nonholonomic) forces. Excellent reviews 
of the history, various forms and geometric descriptions of the reduced systems, 
as well as many relevant examples can be found in |42[ I36L \E\ E| IHl EH , see also 
references therein. 

Apparently, Appel was the first to propose time substitution in order to 
eliminate these extra terms and to transform the reduced systems to a canonical 
(Hamiltonian) form. After that, Chaplygin realized this idea in his reducing 
multiplier theory for nonholonomic systems with two degrees of freedom. 

The key feature in Chaplygins approach is the existence of an invariant 
measure of the reduced system, a rather strong property which puts the system 
close to Hamiltonian ones. For reduced generalized Chaplygin systems emerged 
from classical dynamics, this problem was considered in |2S!- Recently, neces- 
sary and sufficient conditions for the existence of such a measure when the 
Lagrangian of the system is of a pure kinetic energy type are given in |16| . 

On the other hand, numerous attempts to extend Chaplygins reducing mul- 
tiplier theory to systems with more than two degrees of freedom (even having 
an invariant measure) were ineffective, since in this case several overdetermined 
conditions on the metric and constraints are imposed (OOOl). To our knowl- 
edge, until recently there were very few nontrivial examples of multidimensional 
systems that are reducible to a Hamiltonian form exactly by the Chaplygin pro- 
cedure (I^ISniES!). We also quote the results of ^|E!i where some clasical 
nonholonomic problems were reduced to Hamiltonian flows with respect to a 
nonlinear Poisson bracket. 

As an alternative, much effort has gone into the development of the symplec- 
tic and Poisson view of reduced generalized Chaplygin systems. In particular, in 
the case of Abelian symmetries, Stanchenko 44; showed that reduced systems 
can be represented in a Hamilton-like form with respect to an almost symplec- 
tic 2-form 12, which is generally not closed. This observation was extended for 
generic symmetries (see 13 123 El)- In this framework, the Chaplygin multiplier 
is a function / such that the form fH is closed (see |44l 1291 ITC] '). 
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The importance of the existence of an invariant measure for integrabiUty 
of nonholonomic systems was also indicated by Kozlov in '38" 39' , where var- 
ious examples were considered. In ^3 0H1 , Veselov and Vcsclova, inspired by 
classical problems of nonholonomic dynamics, studied nonholonomic systems 
on unimodular Lie groups with right-invariant nonintegrable constraints and 
a left-invariant metric (so called LR systems), and showed that they always 
possess an invariant measure whose density can be effectively calculated. In 
particular, the motion of a rigid body around a fixed point under a nonholo- 
nomic constraint (projection of the angular velocity to the fixed vector in space 
is constant) is described by an integrable LR system (07]). 

Another method of constructing non-Lagrangian (so called L-l-R) systems 
with an invariant measure on Lie groups was proposed in \27\ . The kinetic 
energy of such systems is given by a sum of left- and right-invariant metrics on 
the group. It appears that some L+R systems have natural origins in classical 
nonholonomic mechanics. 

For related problems concerning the integrability of nonholonomic systems 
one can see 03 EH EHl EHl El E3 EH and references therein. Also, the 
existence of an invariant measure for a class of nonholonomic systems with 
symmetries, which include the Chaplygin systems, is recently studied in |52| . 

Contents of the paper. We study several new geometric aspects of non- 
holonomic LR systems on a compact Lie group G. In Section 2 we show that a 
class of such systems can be naturally considered as generalized Chaplygin sys- 
tems on the principle bundle G ^ Q — G/H, where H is a. subgroup of G. Such 
systems are reduced to non-Hamiltonian equations on the cotangent bundle of 
the homogeneous space Q. The latter are described by a Lagrange-d'Alembert 
equation with extra nonholonomic terms which are explicitly found. 

In Section 3 we describe the invariant measure of the original and reduced 
LR systems. If the homogeneous space is two-dimensional, then, by the Chap- 
lygin reducibility theorem, the existence of such a measure leads to a time 
substitution such that our system becomes Hamiltonian. On the other hand, 
we prove that if the reduced system is transformable in this way to a Hamil- 
tonian form for any dimension, then it must have an invariant measure whose 
density is prescribed by the corresponding reducing multiplier. We also show 
that the reduced LR system on Q always possesses an invariant measure; this 
does not necessarily holds for generic Chaplygin systems. 

As a natural example of LR systems. Section 4 describes the classical Veselova 
problem on the motion of a rigid body with a nonholonomic constraint and some 
of its integrable perturbations, as well as its relation to the Neumann system 
and an integrable geodesic flow on the 2-dimensional sphere. 

In Section 5 we consider multidimensional Veselova nonholonomic systems 
on the Lie group SO{n) characterized by various types of constraints and de- 
scribe their invariant measure. The constraints allow a reduction of these sys- 
tems to non-Hamiltonian flows with an invariant measure on the cotangent 
bundle of Stiefel varieties V{r,n). 

In Section 6 we concentrate on the case r = 1, which corresponds to reduced 
flows on the unit sphere S*""^. We show that for a special choice of the inertia 
tensor and after a time substitution, the flow reduces to a completely integrable 
geodesic flow on the sphere. This provides a first example of a nonholonomic 
system with more than two degrees of freedom for which the celebrated Chap- 
lygin reducibility theorem is applicable for any dimension. 
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Also, we prove that, under another time substitution, the multidimensional 
Veselova nonholonomic system on SO{n) reduces to the Neumann system on 

In in last section, for the above integrable case, we explicitly solve the re- 
construction problem: given a trajectory of the reduced geodesic flow on S'"~^, 
to find the corresponding nonholonomic motion on the group SO{n). To per- 
form this, we use the remarkable relations between the Neumann system, the 
geodesic flow on an [n — l)-dimensional ellipsoid, and the evolution of orthog- 
onal frames associated to the geodesies. It appears that the right-invariant 
distribution D C TSO{n) is foliated with invariant tori of generic dimension 
n — 1 and the unreduced LR system is integrable. 



2 Generalized Chaplygin and LR systems on 
Lie groups 

Suppose we are given a nonholonomic Lagrangian system (M, I, D) on the n- 
dimensional configuration space M with (local) coordinates x and Lagrangian 
l{x,x) in the presence of a fc-dimensional distribution D C TM describing 
kinematic constraints: a curve x{t) is said to satisfy the constraints if x(t) G 
for all t. The trajectory of the system x{t) that satisfies the constraints 
is a solution to the Lagrange-d'Alembert equation 

{ dl d dl \ , , 

(----,,j=0, for an ,^D.. (2.1) 

Here (•, •) denotes pairing between dual spaces. 

Now assume that M has a bundle structure tt : M Q with a base manifold 
Q and the map tt is a submersion, that is, TxM = Dx ® Vx for all x. Here Vx is 
the kernel of TxTt and it is called the vertical space at x. Then the distribution 
D can be seen as a collection of horizontal spaces of the Ehresmann connection 
associated with tt : M ^ Q. 

Given a vector Xx G TxM, we have the decomposition Xx — X'^ + XJ^, 
where X^ G Dx, X"^ G Vx- The curvature of the connection is the vertical 
valued two form B on M defined by 

B{Xx,Yx)^~[X^,n]: 

where X and Y are smooth vector fields on M obtained by extending Xx and 
Yx. 

With the help of the Ehresmann connection, the equations of motion can 
be put into the form (see [S]) 

where lc{x,x) = l{x,x^) is the constrained Lagrangian. 

The form of equations H2.2() is very useful in the presence of symmetries of 
the system. Namely, suppose that the configuration space is a principal bundle 
■K : M ^ Q — Mj © with respect to the (left) action of a Lie group (5, and is a 
principal connection (i.e., D is a ©-invariant distribution). Let the Lagrangian 
I also be ©-invariant. Then equations (|2.2|) are ©-invariant and induce a well 
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defined reduced Lagrange-d'Alembert equation on the tangent bundle TQ = 
D/<&. The system {M,l,D) is referred to as a generalized C'haplygin system 
(see EElEl)- 



LR, systems. Now let M be a compact connected Lie group G of dimension 
n with local coordinates g, and g = TjdG its Lie algebra with comutator [, ]. 
Let (•, •) denote either the Adfj-invariant scalar product on g or the bi-invariant 
scalar product on G, and dsj denote the left-invariant metric on G given by 
the nondegenerate inertia operator X : g ^ g in the usual way: 



where uJi — g^^rji, uJ2 = g^^V2- 



Let j/i , . . . , y„ be independent left-invariant vector fields on G generated by 
some basis vectors Yi, . . . ,Yn in the algebra. Following [471 US] , one can define 
an LR system on G as a nonholonomic Lagrangian system (G,l,D) where 
I = ^{g,g) — v{g) is the Lagrangian with a left-invariant kinetic energy and D 
is a right-invariant (generally nonintegrable) distribution on the tangent bundle 
TG. 

The right-invariant distribution is determined by its restriction d to the 
Lie algebra as follows: Dg = d ■ g — g ■ {g^^ • 5 • C TgG, d =const. Let 
f) = span ([)i, . . . , f),„) be the orthogonal complement of with respect to (•, •) 
and f)s =const. Then the right-invariant constraints can be written as 

uieg-'^-d-g, or = (w, • f)^ • g) = 0, s = l,...,m, (2.3) 

where w = g~^ ■ g. 

The LR system (G, /, D) can be described by the Euler-Poincare equations 
(also referred to as the Poincare-Chetayev or Bolzano-Hamel equations) on the 
product X G, 

d ^ 

—lu; = [Iuj,uj] - y{v{g)) + ^Xsg^^ • f)s • , 

9 = ffw, 

where y{v) = {yi{v), . . . , yn{v))'^ is the vector of Lie derivatives with respect to 
the above left-invariant fields j/i, . . . , y„, and are indefinite multipliers which 
can be found by differentiating (|2.3|) . 

These equations define a dynamical system on the whole tangent bundle 
TG, and the right-invariant constraint functions fs in H2.3|l are its generic first 
integrals. Thus, the LR system (G, I, D) itself can be regarded as the restriction 
of the system H2.4|l onto D C TG. (Also, the LR system with non-homogeneous 
right-invariant constraints fs = Cs ^ can be considered as a subsystem of 

For the case v{g) — 0, system (|2.4ll can be reduced to the form 

d ^ 
—Xuj ~ \Xuj, Lj] +y AcJ%, 

dt ^ ' ^ (2.5) 
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where Ts{g) ~ dfs{uj,g)/duj = g ^ • t)s • <?■ This forms a closed system on the 
space {uj,J^i, ■ ■ ■ ,J^s)- 

There is another way of describing LR systems, which is based on the non- 
holonomic version of the Noether theorem (see, e.g., |21 ESj 3.)- Namely, 
as shown in 48 , for v{g) — 0, equations (|2.4() have the conservation law 
^ prg((7 • Tu! ■ g~^) — 0, which can be rewritten as 

^(prg-i.Dg^^) = [Piff-i.os^:^:'^]- (2-6) 

Next, for the case of non-homogeneous constraints fs = Cs, one has 
^(pi'i,(5t^5~^)) = 0, which implies 

^(PV^Dff^) = iP^a-^g^^"^]- 

Combining the above equations, we obtain the momentum equation 

M = [M,uj], (2.7) 

M =prg-i.D-g2:^ + Prg-i-h-s^- (2-8) 

As follows from (|2.8|l . the linear operator sending cj to is nondegenerate, and 
one can express lo in terms of M and the group coordinates g uniquely. Thus 
H2.7|l together with the kinematic equations g — guj represent a closed system 
of differential equations on the space {uJ,g) or on (A^, g), which is equivalent 
to system (|2.4ll . Since on I? C TG we have pi'g-i.o.g = pr^-i.^.^ Xo;, on this 
subvariety the system has the kinetic energy integral ^{M,uj) = ^{Tuj,uj). 

Now let 

t) = span (iTi, . . .,«„_„), (mfc, Ws) = Sks 

and put Wk — g^^ ■ tVk ■ g- Then the above system leads to a closed system of 
differential equations on the space (w, Wfc) or on {M,yVk), 

M = [M,iu], m^[Wk,iu], (2.9) 

n—m 

fe=i 

The distribution D is represented as an invariant subvariety of (|2.9|) given by 
the condition 

n—m n—m 
k=l k=l 

Reduction. Let the linear subspace f) be the Lie algebra of a subgroup 
H G G. Furthermore, we suppose that the potential v(g) is i?-invariant. Then 
the Lagrangian / = ^{g,g) — v{g) and the right-invariant distribution D are 
also invariant with respect to the left iJ-action. (Notice that for m > 1 the 
constraint functions fs themselves may not be 77-invariant.) In this case the LR 
system (G, /, D) can naturally he regarded as a generalized Chaplygin system. 
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Consider the homogeneous space Q = H\G of left cosets {Hg}. The distri- 
bution D can be regarded as a principal connection of the principal bundle: 

H — > G 

i 

Q = H\G 

The Lagrange-d'Alembert equation (|2.2|) is i/-invariant and it reduces to a 
second order equation on Q. In order to write the reduced equations in a simple 
form, we identify q and q* by the -invariant scalar product (•, •), and the 
spaces TQ, T*Q by the normal metric induced by the bi- invariant metric on 
G. Next, consider the momentum maps 

^■.TG^T*G^Q, <P:TQ = T*Q^Q 

of the natural right actions of G on T*G and T*Q, respectively. We have 

(j){X) =g-'-X, Xe TgG 

and the map <& can be considered as a restriction of </) to I?. 

The reduced Lagrangian is by definition the constrained Lagrangian 

lc{9,9) = \{PT^g-^T,g^W9i9)),(l>i9,9)) -v{q) 

considered on the orbit space H\D = T{H\G). It follows that the reduced 
Lagrangian is simply given by 

Liq,q) ^^(miq,q),^q,q)) -V{q), 

where q = '^{g) are local coordinates on Q (which may be redundant) and 
V{q) = v{g). For v = V = 0, this Lagrangian describes a metric which we shall 
denote by dsj jj . 

The reduced system on TQ is defined by the following proposition, which 
appears to be a special case of the general nonholonomic reduction procedure 
described in (3^1 El- 

Proposition 2.1 The reduced Lagrange-d'Alembert equation describing the mo- 
tion of the LR system (G, I, D) takes the form 

(2.10) 

where pr^-ifj^ : g — > g ^fjg is the orthogonal projection and q — Tr{g). 

As a result, l|2.1()|l leads to a system of Lagrange equations on TQ with some 
extra terms. Note that this system always has the energy integral 

E{q,q) ^^(I^q,q),^q,q)) + V{q). 

Proof of Provosition XZ.W First, we need to describe the curvature of the princi- 
pal connection associated to the distribution D. Let Xi,X2 E TgG. Then the 
horizontal and vertical components of Xi have the form 

X!' = 9- Wg-iog (l^iXi), X^ =g- pr^-iijg (j^iXi). 
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Next, if Xi,X2 are right invariant extensions of Xi and X2, then [Xi,X2]g = 
—g ■ [(l){Xi),<j>{X2)]- (Here the first square brackets denote the commutator of 
vector fields, and the second ones represent the commutator in the algebra g.) 
Thus, the curvature is 

S(Xi,X2) = -[Xi',X2'']^=5.pr^-,^^[pr^-,,^0(Xi),PVio,0(X2)]. 

Therefore the right hand side of (|2.2|) is equal to 

(^^'9 ■Prg-it,g[uj,(j){r])]^ = {Iuj,pi-g-n^g[u;,(j){ri)]), uj = g'^ ■ g = <j>{g). 

Combining the above expressions, we arrive at the right hand side of H2.1()|l . 



Reduced momentum equation. Similar to the original LR systems, in 
the absence of potential forces one can describe reduced LR systems on T*Q in 
terms of a momentum equation as well. 

Namely, let us now identify {p} — T*Q and {q} E TqQ by the metric dsj ^, 
i.e., we put p ~ dL{q, q)/dq, and also identify the spaces g — {uj} and g* = {M} 
via relation (|2.8|) . 

Next, introduce the momentum map <i>* : T*Q — > g*, {q,p) A4 hy setting 

^*{q,p)^^q,dL{q,q)/dq), 

where dL{q,q)/dq is considered as an element of TqQ (via identification given 
by the normal metric). 

The map is correctly defined because 

^*iq,p)\p=dLiq,q)/dq = P^g- 1 Dff ^^^(l?, q) =pTg-i^gM. (2.11) 

Indeed, a preimage of dL{q,q)/dq in Dg C TgG, g G Tr^^{q) can be chosen in 
form dlc{g, g)/dg, TT*{g) = q. Therefore, we have 

$(g, dL{q, q)/dq) = pig-i^g (j^ig, dlc{g, g)/dg) 

= Prg-iDg(5 \=Wg-il,g9 [9^9 )5 = PI-g-iDg^^^ 

which establishes the first equality in (|2.1H) . The second equality follows from 

(1231). 

Since the linear subspace g^^dg C g is iJ-invariant, it depends only on 
q G Q. Thus, the system (|2.9|) represented in terms of w can be regarded as 
a flow on the quotient manifold H\TG = Q x g obtained from TG = G x g 
by factorization by H. The same system represented in terms of A4 leads to a 
system on Q x g*. 

Relations between the above manifolds are described by the commutative 
diagram below, where the vertical arrows denote the corresponding inclusions 
and $* are the extensions of the momentum maps $, $* respectively. 

H\ (2.8) 
TG ^ Gxg U g X g ^ U Q x q* 



D = G X D > I Q > 1 Q 
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For a fixed q, the map $* establishes a bijection between the subspace g^^tig C 
0* and the cotangent space T*Q. 

Now, applying (|2.6(l and (|2.1HI . we arrive at the reduced momentum equa- 
tion 

j^^*{q,p) = [^*{q,p),^{q,q)], (2.12) 

where q — q{q,p) is determined from p = dL{q,q)/dq. This leads to a system 
of equations on T*Q, which are equivalent to the Lagrange equations on TQ 
obtained from H2.10|l . 

As a consequence of the momentum equation H2.12II . we also obtain the 
following result. 

Proposition 2.2 In the absence of potential forces the reduced LR system on 
T*Q always has a set of first integrals A = {f o / e M[g]'^}, where R[3]^ 
is the algebra of Ada invariants on g . 

The number of independent functions in A is equal to the number of inde- 
pendent G-invariant functions on T*Q, that is to dimpr[,(ann(^)), for a generic 
^ G c) (see [^). Here ann(^) = {jy g g, rj\ — 0}. \i Q = H\G is a symmetric 
space, this number is equal to the rank of Q. 



3 Invariant measure and time rescaling 

One of the remarkable properties of LR systems is the existence of an invariant 
measure, which puts them rather close to Hamiltonian systems. 

Theorem 3.1 (gTlllHj). The LR system ifF3)) on the space {uj,Ti, . . . ,Ts) 
possesses an invariant measure with density 

fi = y/dct(X-i|g-i|,g) = y^det{J^s,T-^J^i), s,l = l,...,m, (3.1) 

where 2^^\g-itjg is the restriction of the inverse inertia tensor to the linear 
space g~^i)g C g. 

The alternative description of LR systems given by the momentum equation 
H2.7|) leads to another expression for invariant measure. 

Theorem 3.2 The LR system defined by the momentum equation ^2.9\) has 
the invariant measure 

fidioA dWi A • • • A dWn-m = m"^ dM A dWi A • • • A dWn-rn , (3.2) 
= ^det(X|<,-ia,) = y^dct(W„XW,), (3.3) 
i,j = l,...,n-TO, 

where X\g-i^g is now the restriction of the inertia tensor to the linear space 
g-'^dg c g. 

Expressions (|3.1|l and l|3.3|l involve complimentary basis vectors in g^^29- 
this sense the densities fi and fi given by the above theorems are dual. 

Proof of Theorem \3./A First note that the systems (|2.4|) and 1)2. 8|) can be ex- 
tended to one and the same system on the space (w, . . . , Tm, Wi, ■ • ■ , Wn-m) 



dM 
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by adding evolution equations for and Ti respectively. The resulting sys- 
tem has an invariant measure whose density can differ from those of the original 
systems only by constant factors. Hence the functions ^ in %\.\\ and jl in (|3.2|l 
can be different only by a constant multiplier. 

Next, note that in an appropriate g-dependent orthogonal basis in the alge- 
bra g, the Jacobian matrix dM/doj has the following block structure 

dM ^ 0\ /O \ (I\g-^^g S \ 

where In-m i Im are unit matrices of dimension (n — m) x (n — m) and m x m 
respectively, and iS is some {n — m) x m-matrix. In the same basis one has 

dM . ^ 0\ /O \ 1 ^ (In-m \ 

with some m x [n — TO)-matrix hi. Comparing the right hand sides of these two 
expressions with we obtain the following chain: 

= det(X-i)det(X|5-iag). 

(3.4) 

Hence, we can choose the density jl in the form (|3.3f) . 
Finally, taking into account the relation 

-1 

dM A dWi A • • • A dWn-m 

and using we come to the equality in H3.2|l . The theorem is proved. 

As shown in ^S] , Theorem l3 . 1 l implies that the original nonholonomic system 
H2.5|l on the left trivialization g x G of TG has the invariant measure fi{g) duj/\dg. 

Reduced invariant measure. Now we proceed to reduced LR systems. 
As a natural consequence of the above theorems, we have 

Theorem 3.3 The reduced LR system f2.10\) (or, after the Legendre transfor- 
mation, the system 1^3.1^) on T*{H\G)) possesses an invariant measure. 

Note that the reduction of a generic Chaplygin system may not have this 
property (see P^l. 

The proof of Theorem 13.31 consists of two steps. First, it is seen that the 
restriction of the LR system (|2.5ll onto the distribution D C TG has an in- 
variant measure. Indeed, the volume form on the tangent bundle admits the 
decomposition 

dujAdg^ 9{g) d/i A • • • A d/™ A 77 , (3.5) 

where fs{g,g) are the constraint functions in 12.3|l . 9{g) is a function, and 77 is 
a volume form on D. Since the 1-forms dfs are independent on TG, 9{g) does 
not vanish on G. 

Let be the Lie derivative with respect to the nonholonomic flow H2.5|l . 
Since the functions fsig,g) sue its generic first integrals, we have C^d fs = 
d{fs) — 0, s — 1, . . . ,m. As a result, from the condition £,(/i dio A dg) — and 



det(X ^|g-l|,g) 



dM 



duj 



det(X-i) 



dM 



duj 



dw A dWi A • • • A dWn 



dM 
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1)3.5(1 we obtain d/i A • • • A d/m C.* (m^ = 0. Hence, the restriction of the flow 
onto D has the invariant measure fJ.ig)9{g) U. 

Notice that one can always choose U to be i/-invariant. In this case, since 
the form duj A dg is G-invariant, whereas the wedge product d/i A • • • A dfm 
and fi(g) are _ff-invariant, the density fJ'{g)d{g) of the restricted measure is also 
7J- invariant and goes down to Q. 

The second step is based on the following general lemma. (Although it is 
quite natural, we could not find it in the literature.) 

Lemma 3.4 Suppose a compact group © acts freely on a manifold N with local 
coordinates z, and there is a (5-invariant dynamical system z = Z(z) on N. 
If this system has an invariant measure (which is not necessarily (5-invariant), 
then the reduced system on the quotient manifold N/& also has an invariant 
measure. 

Now, identifying the group and the manifold N with H and D respec- 
tively, we arrive at Theorem 13.31 

Proof of Lemma 13.41 The manifold N can be locally represented as a direct 
product M'"'{x} X (5, where x is a local coordinate system on N/<3, so that the 
©-action and the dynamical system take the form 

a - {x,g) = {x,ag), a e (S, 
and x^X{x), g = Y = g-^{x), ^(a;) e g = T/d©, 

respectively. 

Let Q be an invariant measure of the original system on N, fi he a bi- 
invariant volume form on 35, a be a volume form on iV/(5 and ax be its local 
representation in ^-coordinates. Then the invariant measure on N locally has 
the form 8 = /(x, g) ^ A Ux- Thus 

l^zf{x; g) A cT, = Z{f) fiAax + f- (Czfi) A ct, + A {Cz(Jx) = 0, (3.6) 

where Cz is the Lie derivative with respect to the flow Z. Since dax = dx<Jx — 
and d/i = dg/i — 0, we have 

Cz<Jx ^ {do iz)ax + {iz ° d)(Tx = d[iz(Jx) = dx{ix(^x) = C^x^x, (3.7) 



Lz[i. = {do iz)n + {iz o d.)fj, = d{izfJ-) = d{iYf^) = {dx + dg){iYfJ.)- (3.8) 

For a fixed x, Y — Y{x) is a left-invariant vector field on 0, whereas the 
corresponding flow on © is right-invariant. Since /i is bi-invariant, we have 
CyfJi = dg{iYlJ-) — 0. Also, it is obvious that dx{iYfi) Act = 0. Therefore, taking 
into account (|3.6II3.8|I . we get 

Z{f)^lAax-^f^l/\{Cx<Jx) = (3.9) 

Now we introduce the "averaged" density f{x) = f{x,g)fi, which, as we 
shall see below, has the following property 

L ^^^^^ ^^{L " ^^'^^^ 
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Then, by integration of H3.9|l . we obtain X{f)(Jx + fCx^x = 0. As a result, the 
reduced system preserves the volume form f{x)ax- 

We stress that the above procedure does not depend on the choice of the 
local coordinates on N/(3. Indeed, let y — y{x) be another coordinate system. 
Then 

6 = Hy,9)f^^crv = h{y{x),g)^i Add CTj, = J{x,g)^iA a^, 

and after integration we have f(x)ax = h{y)ay. 

It remains to prove (101 . We have Z{f) = X{f) + Y{f) and X{f)^^ = 
X (Jig //i) . Therefore the relation H3.10|l is equivalent to 

'rav^o. (3.11) 

To check the latter relation, we fix x. Then Cyiflj) = y{f)l-t + f^YfJ' and, 
on the other hand, Cyiflj) = dg{iY{ffij)- Since £y/x = 0, we get Y{f)n, = 
dgiivif fJ')), and l|3.11|l follows from the Stokes theorem. The lemma is proved. 

Chaplygin reducing multiplier. Here we continue with the reduced 
LR systems. However, all considerations hold for an arbitrary generalized 
Chaplygin system with the Lagrangian of the natural mechanical type. Let 
5i, . . . , gfc be some local coordinates on the homogeneous space Q and pi, . . . 
Pi = dL/dqi be canonically conjugated momenta which provide coordinates on 
the cotangent bundle T*Q. Let gij denote metric tensor of c?Sj ^ and g^^ give 
the dual metric on T*Q. 

The reduced Lagrangian is L{q, q) — ^ Tldijiiqi ~ ^{q)- We also introduce 
the Hamiltonian function H{q,p) = ^Y^9^''PiPj + ^{q) (the usual Legendre 
transformation of L). Then (|2.10() can be rewritten as a first-order dynamical 
system on T*Q: 

qi = —, pi = -- hnj(q,p), i = l,...,fc. (3.12) 

OPi oqi 

The functions H,; are quadratic in momenta and can be regarded as non- 
Hamiltonian perturbations of the equations of motion of a particle on Q. 

Now consider time substitution dr — Af{q)dt, where Af{q) is a differentiable 
nonvanishing function on Q, and denote q' = dq/dr. Then we have the following 
commutative diagram: 

TQ{q,q} ^^^^ TQ{q,q'} 
T*Q{q,p} T*Q{q,p}. 

The Lagrangian and Hamiltonian functions in the coordinates {q, g'} and 
take the form 

L*{q. i')-\Y. ^'aMj - y{q): H*{q,p) = ^Y. -j^a^'p^Po + v^q)- 

There is a remarkable relation between the existence of an invariant measure 
of the reduced system H3.12|l and its reducibility to a Hamiltonian form. 
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Theorem 3.5 1). Suppose that after the time substitution dr = M{q)dt the 
equations J^SA'^) become Hamiltonian, 

Then the original system J^S.li^) has the invariant measure 
N{qf-^ dpi A • • • A dpfe A dgi A • • • A dqfe = N{qf~^ 

where fi = "^dpi A dqi is the standard symplectic form on T*Q. 

2). For k — 2, the above statement can also be inverted: the existence of the 
invariant measure with the density JV{q) implies that in the new time 
dr = J\f{q)dt, the system I^S.l^) gets the Hamiltonian form ]S.13\) . 

In nonholonomic mechanics the factor A/" is known as the reducing multiplier, 
item 2) of this theorem is referred to as Chaplygin's reducibility theorem (see 
[nilTHlIini or section III. 12 in gg). Item 1) of the theorem was impUcitly 
formulated in ^ Ej . 

Remark 3.1 Notice that for k > 2, the multiplier Af{q) and the density of the 
invariant measure of system H3.12|l do not coincide. 

On the other hand, the paper 'TT gives examples of reducibility of nonholo- 
nomic systems to Hamiltonian ones with respect to nonlinear Poisson brackets. 
In some of these examples, even for k — 2 the reducing factor and the density 
of invariant measure are different. 

Proof of item 1) of Theorem \'6.5\ For simplicity we shall use the vector notation 
p = {pi, . . . ,Pk), q — {q^, ■ ■ ■ , q'^), etc. Let G be the matrix (g*-'). Then q = Gp, 
H=l{Gp,p), H* ^^{Gp,p). 

The equations (|3.13l) in the original time t take the form 

q = JVVpH*iq,p), f>=-NV,H*[q,p). (3.14) 

Equations H3.12|l have an invariant measure with density / if 

(V„ /VpiJ) + (Vp, f{-V,H + n)) = 0. (3.15) 

For / which depend only on q-coordinates, we have (Vp, n) + (Vg In /, Gp) = 
0, or equivalently 

d{\nf)+a = d{\nf) + {A,dq)^Q, (3.16) 

where (Vp,n) — {A,q) — a{q). In particular, the one-form a is closed. 

We shall prove that the function f{q) = M^~'^{q) satisfies equations (|3.16|l . 
Since p = Mp we have Mp + Np = p. Therefore, using equations l|3.14(l we 
obtain 

p = M-'p - N{q)p = -V,H*{q,p) - {V,N, Gp)p. (3.17) 

Also, one can easily see that '\/qH*{q,p) = W qH {q, p) - JV^^ {Gp, pjW qAf . Thus, 
comparing H3.12|l and H3.17|l we get 

Il{q,p)^^-\Gp,p)Vq^~Af-\VN,Gp)p (3.18) 
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Using l|3.18|l we see 



(Vp, n) = -1 (2(V,AA, Gp) - fc(V,AA, Gp) - (V,AA, Gp)) - ^(V,AA, Gp). 

Hence a = ~d\n{J\f''-^). 

As mentioned above, item 2) of the theorem is just a reformulation of Chap- 
lygins reducibihty theorem in 

Clearly, the density of an invariant measure of a generic dynamical system 
depends on the choice of local coordinates on the phase space. However, in the 
case of a system on the cotangent bundle T*Q the density is invariant with 
respect to changes of coordinates on Q, since the symplectic form f2 and the 
measure itself are invariant with respect to contact transformations. 

Remark 3.2 The paper ^H* (see also |TQ) contains a nontrivial observation 
about the density of the invariant measure, which in our terms reads as follows. 
Let a function f{q,p) be a solution of H3.15|l in case of the absence of the 
potential {V{q) = 0). Then one can check that the function fo{q) = /(g, 0) 
also satisfies condition H3.15|l . i.e., it is a solution of H3.16() . In other words, if 
the reduced system (I3.12f) has an invariant measure for V ^ 0, one can take 
this measure to be of the form f{q)f2'^. Then, since (|3.1()|) does not depend on 
the potential, the reduced system H3.12() has the same invariant measure in the 
presence of the potential field V{q) as well. 

4 Veselova system on T SO{3), the Neumann sys- 
tem and a geodesic flow on S'^ 

The most descriptive illustration of an LR system is the Veselova problem on the 
motion of a rigid body about a fixed point under the action of the nonholonomic 
constraint 

(f^,7) = 0. (4.1) 

Here fl is the vector of the angular velocity in the body frame, 7 is a unit 
vector which is fixed in a space frame, and ( , ) denotes the scalar product in 
K.^ 021 ■ Geometrically this means that the projection of the angular velocity 
of the body to a fixed vector must zero. 

This setting should not be confused with the nonholonomic Suslov problem, 
when the analogous constraint is defined by a vector fixed in the body frame 

(|sniEHlE3)- 

The equations of motion in the moving frame in the presence of a potential 
field V = 1^(7) have the form 

dV 

in = If7 X f7 + 7 X — + A7, 
07 

7 = 7x1], (4.2) 

where J is the inertia tensor of the rigid body, x denotes the vector product in 
K^, and A is a Lagrange multiplier chosen such that 0(t) satisfies the constraint 

(J ^7,7) 



14 



The Veselova system (I4.1|) . (|4.2() is an LR system on the Lie group 50(3), 
which is the configuration space of the rigid body motion. After identification 
of the Lie algebras (R'^, x) and (so(3), [•, •]), the operator X induces the left- 
invariant metric ds^. The angular velocity vector 17 correspond to g^^g, the 
velocity in the left trivialization TSO{i) = SO{3) x so(3), and the Lagrangian 
function equals ^ (17,117) — 1^(7). The fixed vector in the space corresponds to 
the right-invariant vector field "fg — g ■ {g^^ • t) • g) G TgSO{3), () G so(3), and 
the nonholonomic constraint H4.1|l has the form {g~^ • [} • 17, 17) = 0. 

Equations (|4.2I) . H4.3|l also define a dynamical system on the space {17, 7} = 
so(3) X R^, and the constraint function (17,7) appears as its first integral. As 
noticed in |47| . this system has an invariant measure with density ^ {T^^j, 7). 
Apart from the above constraint, it always has the geometric integral (7,7). 
According to 26 , for ^(7) — there are two other independent integrals 

i(17,X17)-(17,7)(I17,7), i(X17-(X17, 7)7-^(17,7)7)2, (4.4) 

the first expression being an analog of the so called Jacobi-Painleve integral 
which replaces the energy integral in some systems with nonstationary con- 
straints. 

On the constraint subvariety (|4.1|l . these functions reduce to the energy 
integral Fi = ^(X17,17) and an additional integral 

^^2 = i(/17,/17)-i(/17,7)2, 

which was found in |47| . 

As a result, by the Euler-Jacobi theorem (see e.g., |2|), the above system is 
solvable by quadratures on the whole space so(3) x R'^. Note that analogous 
integrable LR systems on the group S'i(2,R) and the Heisenberg group are 
studied in 

As shown in in the case of the absence of the potential the Veselova 
system 1)4.211 . (|4.1|l can be explicitly integrated by relating it to the classical 
Neumann system. 

Theorem 4.1 (@H1)- 7(*) o. solution of equations \4-''^ , \4-l^ with 

^(7) — and with energy constant Fi — h. Then after time reparameteri- 
zation 

/2/idetI-i , 

the unit vector q = j is a solution of the Neumann system on the unit sphere 
S'^ = {q GM.^ \ qf + qj + q1 = 1} with the potential U{q) = i(Ig, q), 



drf 

corresponding to the zero value of the integral 



^q = -Xq + \q, (4.5) 



I ( -^q X q] , -^q x g ) - detX (I^^g, q). (4.6) 



dri J dri 

We emphasize that for (17,7) 7^ 0; Theorem 14 . 1 1 does not hold, and in this 
case the procedure for integrating equations (|4.2() . H4.3|l was indicated in |26| . 
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Reduction. The above relation between the LR system and the Neumann 
system via the time reparameterization appears to be quite natural in view of 
the fact that the Veselova system is a Chaplygin system on the S'0(2)-bundle 

50(2) — > 50(3) 

I 

52 = 50(2)\50(3) 

where SO{2) is the subgroup generated by rotation about the vector 7. Indeed, 
the Lagrangian and the nonholonomic constraint H4.1|l are invariant with respect 
to such rotations. Hence, the Veselova system can be reduced to the (co)tangent 
bundle of = {q e \ qj + ql + ql ^ !}■ 

The momentum map $ : TS'^ so(3) = is simply given by ^{q,q) = 
q X q, hence the reduced Lagrangian is L{q,q) — ^ {l{q x q),qx q) — V{q). 
Note that the reduced potential is given by the same function V, regarded as 
a function of q instead of 7. 

Next, in view of the relation 

prff-ii,ff [^(9. 9), C)] = (9, (9 X 9) X (C X q))q = 9 x 

where ^ = (^11^2,^3)"^ is any tangent vector of 5^ at the point q, the reduced 
Lagrange-d'Alembert equation H2.1()|l takes the form 

^-^^,ej=*(<Z,.,a ^-HiqxqUxO. 

Now the reduced LR system on T*S^ can explicitly be written in terms of 
local coordinates 91,92 on and the corresponding momenta pi = dL/dqi, 
P2 = dL/dq2, 

dL d d^ , , 

dqk dt de,k 
where L, are obtained from L{q,q), ^'((7,(?,0 by the substitutions 

qiqi + 9292 . <7iCi + ^26 

93 = , „ =^1 43 



- 9i - 92 y/'^-ii- 92 

A direct (but tedious) calculation shows that the reduced system (|4.7|) has 
an invariant measure with density M{q) = 1/ ^ {q,T^^q). (As was mentioned 
above, the latter does not depend on the choice of local coordinates on 5^). 

Since the reduced system is two-dimensional, Chaplygin's reducibility the- 
orem (item 2 of Theorem 13. 5f) says that in the new time dr = J\fdt and new 
momenta pk = Afpk, k = 1,2, equations 14.7|l transform to a Haniiltonian sys- 
tem. Equivalently, the latter is described by the following Lagrangian obtained 
from L{q,q), 

L*iq, q') = ^^^Z^ X q), q' X q)) - V{q), q' = ^. (4.8) 

For = 0, this is a Lagrangian of a geodesic flow on S^. 
Theorem 4.2 The geodesic flow on 5^ with the metric 

{q,X-\)-^ds]^0, ds]^o = detX [{dq,I-^dq){I-\, q) - dqf] 
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obtained from is completely integrable. It has an additional quadratic 

integral 

which corresponds to the integral F2 of the LR system \4--'/J^ , \4-l\l - 

This theorem, as well as our observations on the reducibility of the Veselova 
system to Hamiltonian form, is a part of a general integrability theorem for 
a multi-dimensional Veselova system on the group SO{n), which we discuss 
in detail in Section 6. In the rest of this section we only quote some specific 
properties of the 3-dimensional case. 



The Veselova system on 50(3) with integrable potentials. Clas- 
sical integrable cases of the rigid body motion without constraints were already 
used to produce integrable geodesic flows on the sphere (see, e.g., [01). Namely, 
consider the Euler-Poisson equations of the motion of the rigid body with tensor 
of inertia J7 and axisymmetric potential V(7) 



dV 

jn = jnxn + -f X —, 7 = 7x17, 
07 



(4.9) 



which always have first integrals 



(7,7) = !, Z2-(^r!,7), 



f,^-{jn,n) + vi-f). 



In the Euler case (V(7) = 0) there is an additional integral /2 = ^{Jil, Jil), 
and under the condition i2 — and the substitution 7 = 5, equations (|4.9|l 
define a geodesic fiow on the sphere 5^ with the metric 

rfsj.p = (g^Jg) ^^^"^'^'^'^^' 

There is an interesting duality between integrable potentials and additional 
first integrals of the Euler-Poisson equations and of the Veselova system. 

Lemma 4.3 The Veselova system i4.1\ ), i4.S^ with the potential V{j) has an 
additional integral of the form F — F2 + F{'y) and therefore is integrable by 
the Euler-Jacobi theorem if and only if the Euler-Poisson equations \4-^ with 
inertia tensor J = T^^ and the potential V — F^j) are integrable for i2 — 
due to the presence of the extra integral f2 + V{'j). 

Proof. Indeed, the necessary and sufficient condition for equations (|4.9|l with 
12 = to have the integral ^{J'il, J^O.) + V{-f) has the form 



7 X 



dV 



dV 



X 7 



X 7 = . 



On the other hand, the system has the integral \ (JfJ, 1VL)-\ {IVL, 7)^ 

i^(7) if and only if 



7 X 



dF 

97 



d-f 



X 7 



X 7 = . 
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Since we set J ^ = X, ^(7) = "^'(7), both conditions arc equivalent, which 
proves the lemma. 

Some integrable polynomial potentials for the Euler-Poisson equations are 
given in . In a similar way, one can construct integrable polynomial potentials 
(or Laurent polynomial potentials, such as given in [52]) for the Veselova system. 
For example, the following proposition holds. 

Proposition 4.4 Let I — diag {Ii, 12,13). The Veselova system 
with potential 

V{j) = a, ((X^,7) - (J7,7)') +«2(X7,7) + ^ + ^ + ^, (4.10) 

7i 72 73 

ai, . . . , as being arbitrary constants, is solvable by quadratures. The additional 
integral is: 

F = i(Xr!,l])- i(Xl],7)2+aidetX(X7,7)(X-i7,7)-a2detX(X-S,7) 

+^3 (l2^+h^)+ai (h^+Ii^) +«5 (h^+h^) . 
\ it it J \ 72 72 / V 73 73/ 

Note that the integrability of the Veselova system with the Clebsch potential 
a (17,7) was already shown in |T7ll^ . 

5 Nonholonomic LR systems on SO{ri) and their 
reductions to Stiefel varieties 

Now we proceed to a generalization of the Veselova system, which describes the 
motion of an n-dimensional rigid body with a fixed point, that is, the motion 
on the Lie group SO{n), with certain right-invariant nonholonomic constraints. 

For a path g(t) G SO{n), the angular velocity of the body is given by the 
left-trivialization lolt) = g^^ ■ g{t) G so(n). The matrix g G SO{n) maps a 
coordinate system fixed in the body to a coordinate system fixed in the space. 
Therefore, if ei = (en, . . . , ei„)'^, . . . , e„ = (e„i, . . . , e„„)^ is an orthogonal 
frame of unit vectors fixed in the space and regarded in the moving frame, we 
have 

El ^ g ■ ei, En = g- Cn, 

where = (1, 0, . . . , 0)^, . . . , = (0, . . . , 0, 1)^. From the conditions = 
Ei = g-Ci+g-ei, we find that the vectors ei, . . . , e,i satisfy the Poisson equations 

ei = -Lue.i, i = 1,. . . ,n. (5.1) 

Below we use the convention a; Ay = x®y—y(>^x — x-y"^ —y-x^ . Also now (•, •) 
denotes the Killing metric on so{n), {X,Y) — —^tr{XY), X,Y € so{n). The 
left-invariant metric on SO{n) is given by a non-degenerate inertia operator 
X : so{n) —f so{n) and the Lagrangian of the free motion of the body is 
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Right-invariant constraints on SO{n). What form may have a multi- 
dimensional analog of the classical constraint (|4.1|l ? To answer this question, 
we first note that instead of rotations about an axis in the classical mechanics, 
in the n-dimensional case there are infinitesimal rotations in two-dimensional 
planes spanned by the basis vectors ei,ej, i,j = l,...,n. Suppose, without 
loss of generality, that in the three-dimensional case 7 = ei. Then condition 
H4.1|) can be redefined as follows: only infinitesimal rotations in the planes 
span(ei,e2) and span(ei,e3) are allowed. Hence, it is natural to define its n- 
dimensional analog as follows: only infinitesimal rotations in the fixed 2-planes 
spanned by (ei, 62), . . . , (ei, e„) (i.e., in the planes containing the vector ei) are 
allowed. This implies the constraints 



(w, Ci A Cj) — 0, 2 < i < j < n. 



(5.2) 



Following 28 , one can relax these constraints by assuming that the angular 
velocity matrix in the space has the following structure 



uj = gujg 



( 



LO\r 



o 



-LOr 



where O is the zero (n — r) x (n — r) matrix. 

Equivalently, consider the right-invariant distribution on TSO{n) whose 
restriction to the algebra so{n) is given by 



= span{£^j /\ Ek, /c = 1, . . . , r, j = 1, 



A, 



where Ei A Ej form the basis in soin). Since Ci A ej = g ^ ■ Ei A Ej ■ g, we find 
that the constraints are defined by relations 



ui Vr — g ^ ■ d ■ g — spanjei A e^, , 



A e,;, I < i < n}, 



that is (w, Cp A Gg) = 0, r < p < q < n. 



(5.3) 



The LR system on the right-invariant distribution Z),. C T SO{n) can be 
described by the Euler-Poincare equations (|2.4I) on the space so[n) x SO{n) 
with indefinite multipliers \pq, 



i + uje,. = 



E 

r<p<q<.7i 



^pq ^p ^ ^q: 



i^i,. 



(5.4) 



Here the components of the vectors ei, . . . , e„ play the role of redundant co- 
ordinates on SO{n). For n = 3,r = 1, after identification of the Lie algebras 
(]R.^,x) and (so(3), [•, •]) and setting = Sijk^k, £1 = 7, this becomes the 
Veselova system (|4.2f) with V — 0. 

By analogy with (|4.2|l . we will call (S'O(n), I, Dr) a multidimensional Veselova 
system. 

Differentiating (|5.3|) , from (|5.4I) one can obtain a system of linear equations 
for the determination of the multipliers in terms of the components of ti;, u>, and 
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ej. Thus, (|5.4|l contains a closed system of differential equations on the space 
(wy , Cr+i, . . . , e„). The latter system has first integrals 

(w, Cp A Cq) = Wpq, Wpq = COUSt, T < p < q < U (5-5) 

and our LR system on Dr C T SO{n) is the restriction of H5.4I) onto the level 
variety Wpq = 0. 

As follows from Theorem 13.11 the system (|5.4() has an invariant measure 
with density 



= Vdct = \J\{ep A eg,X-i(e5 A e;))|, 

r < p < q < n, r < s < I < n, 

where l.T>r C so{n) is the orthogonal complement of Vr with respect to the 
metric (•, •) and I~^\±-Dr is the restriction of the inertia tensor to _L2?r- 

In the case of the Veselova system on TSO(3) after identifying ei = 62 x 63 
with 7 S R"^ the above expression reduces to the known form ^ (7, 

In practice, for a large dimension n and small r, the number of constraints 
H5.3(l is large, which leads to rather tedious expressions for the explicit form of 
the system and the density of its invariant measure. In this case one can make 
use of the alternative momentum description (system H2.7|l ). Namely, in view 
of the matrix representation 

\fXeso*{n), prj,^{x) = TX + XT -rxr, 

r = ei (g) ei + • ■ • + Cr (8) , 
system 1)2. 9|l takes the following form: 

M = [M,uj], (5.6) 
M = prp^(Xtj) + pr^p^ oj 

= uj + {Iuj- uj)T + T{Iuj - - r{Iuj - uj)T . (5.7) 

The map ui ^ A4 given by H5.7|l is nondegenerate. As a result, equations H5.6|l . 
H5.7|l together with the Poisson equations H5.1|l . which are equivalent to 

f-[r,H, (5.8) 

represent a closed system of differential equations on the space {uj, ei, . . . , e^) 
or on the space {A4, ei, . . . , e^). 

In the classical case of n ~ 3,r — 1, passing to the vector variables 7, fi, m, 
where Aiij = SijkVnk, we obtain 

m^m-{mr/h + {n,-fh, and m^m - ^^J~^^'^\ 

(7,1-17) 

Substituing this into (|5.t)|) yields explicit equations in terms of fl, 7, which again 
describe the Vcsclova LR system (I4.2|) . (|4.3|) with V" = 0. 

As follows from the structure of the system (|5.6() , (|5.8() , it possesses a family 
of integrals given by nonzero coefficients of the following polynomial in A 

tr (TW + A^)^ keN. 
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In addition, it has the invariant variety defined by the condition 

M A ei A • • • A = w A ei A • • • A = 0, (5.9) 

where M, lu are considered as 2-forms and Cfe as 1-forms in the same EucUdean 
space M". This gives a set of scalar conditions on the components of or w, 
which describes the Unear subspace Dr — g^^ ■ X) ■ g C soin). Hence, among 
conditions (|5.9|l only (n — r)(n — r — l)/2 are independent. 

Next, according to Theorem 13. 21 the LR system (|5.6|) . (I5.7|l . (|5.8|l possesses 
an invariant measure with dual density 

Q ^ jlduj A dei A ■ ■ ■ hdcr ^ jl^^ dM A dei A • • • A de^ 

ft = ^det(X|cJ = y^|(ei A ep,X(ej A e^))], (5.10) 

\ < p < q < 1 < i < j < n, 

where Xju^ is the restriction of the inertia tensor to the subspace Vr C so{n). 

Remark 5.1 Since I?,, is invariant under the action of SO{n — r) on the linear 
space spanned by the vectors Cr+i, . . . , e„, expression (|5.1U|) . in fact, does not 
depend explicitly on the components of these vectors. Moreover, Vr is also 
invariant under the S'0(r)-action on the space span (ei, . . . , er). As a result, 
the above density depends explicitly only on coordinates on the Grassmannian 
G{r,n) = SO{n)/{SO{n — r) x SO{r)), that is, on the Plucker coordinates of 
the r-form ei A • • • A e^, which are invariants of the above actions. A simplified 
expression for the density depends on the choice of T. 



The special inertia tensor. It appears that for some special inertia ten- 
sors, the density (|5.1U|) takes an especially simple form, which we shall use 
in the sequel. Suppose that the operator 2 is defined by a diagonal matrix 
A — diag{Ai, . . . , An) in the following way 



liE^ A Ej 



dctA ^ 



(5.11) 



Notice that for n = 3 this corresponds to the well known three-dimensional 
vector formula T{x x y) = ^^j: ^ Ax x Ay, A — and thus in this case defines 
a generic inertia tensor. 



Theorem 5.1 Under the above choice oJT, 



det(J|i,J = Vn,r = (dct^)'' 



^A,, •••A:,(ei A--- Ae^)f 



(5.12) 



where p is an integer constant, the summation is over all r-tuples 

I = {1 < ii < ■ ■ ■ < ir < n}, and (ei A • • • A e^)/ are the Plucker coordinates of 

the r-form ei A • • • A e^. 

Proof. It is more convenient to calculate first the dual determinant 
which can be represented in the form 



detA{ep,A ^es){eq,A ^ei) - det A{ep, A ^e/)(e,,^ ^e^) 
r < p < q < n, r < s < I < n. 



(5.13) 
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Since we deal with purely algebraic expressions, in this proof we can regard 
er+i,...,e„ as vectors in the complex space C". Next, since the action of 
SO{n — r) on the linear space A C C" spanned by Cr+i, does not change 

_LI?r C A'^C", the above determinant must depend only on the Pliicker coordi- 
nates 



{er+i A • • • A e„)j, J = {ji 



Jn-r}, 1 < il < ■ • • < jn-r < n. 



In view of dimension and the structure of the determinant (|5.13|) , it is a homo- 
geneous polynomial in the components of Cp of degree 

4 • dim SO{n — r) = 2{n — r){n — r — 1). 

Hence, it is a homogeneous polynomial of degree 2(n — r— l)in(er+iA---Ae„)j. 

Suppose that the linear space A is tangent to a (possibly imaginary) cone 
JC — {{X, A~-^X) = 0} C C" and, without loss of generality, assume that e„ is 
directed along the tangent line Afl/C. Then (e„, A~^ep) = for p = r-l-1, . . . , n, 
and in this case the last n — r — 1 rows and columns of the determinant H5.13II . 
and therefore the determinant itself, vanish. 

On the other hand, the condition for A to be tangent to /C has the form 



det(A-i|A) = 



{en,A ^e^+i) 



r+l 



A- 



(^n; A C71) 



where ^~^|a is the restriction of the quadratic form A onto A. Expanding the 
latter determinant we see that it equals A~^ ■ ■ ■ A~^(e,-+i A • • • A e„)j, thus 
it is a quadratic polynomial in the above Pliicker coordinates. 

Combining our considerations, we see that when A is tangent to /C, the 
matrix A~^|a ha-s corank 1, whereas the matrix 2~^\±xir has corank {n — r — 
1). Hence, the determinant (|5.13() is divisible by the [n — r — l)-th power of 
det(A~^|A), a homogeneous polynomial of degree 2{n — r—l) in the coordinates 
{cr+i A • • • A e„)j. Thus the quotient of the determinant and the polynomial 
has zero degree in these coordinates. Since the quotient cannot have poles, it 
is a constant. An additional study of H5.13|l shows that this constant must be 
a positive power of det A. As a result. 



det(X" 



.) = {AeiAY 



-1 



. ,7 



(e 



r+l 



Ae„)^ 



Pi G N. 



Now, in order to obtain the desired expression H5.12|l . we use the rela- 
tions (|3.4|) with g'^X)g = Vr, g~^i)g = l-Vr, and (e,.+i A • • • A e„)j = (ei A 
• • • A er)j, where / and J are complimentary multi-indices in the sense that 
{zi, . . . ,ir, ji, . . . , jn-r} IS a permutation of {1, . . . , n}. This, together with the 
fact that detX is a power of det A, proves the theorem. 

From Theorems 13.11 15.11 we get 

Corollary 5.2 Under the condition 1^5. the LR system H5.6() - (|5.8(l has an 

invariant measure 



^A,, ■■■Ai^{ei A ■■■Aer)j 



-(ri-r--l)/2 



dM Adei A ■ ■ ■ Ader . (5.14) 
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In the particular case of r = 1, the density fl in ^5.10\) is proportional to 
(ei,v4ei)("-2)/2. 

As follows from (|5.12|) or (|5.14|l , in the opposite extreme case r = n — 1 (no 
constraints) Vn.r and fl are just constants, as expected. 

Reduction to Stiefel varieties. Now we notice that in the case of the 
constraints (|5.3|l , the orthogonal complement f) of i) is a Lie algebra, namely 

f) = span{_Ep A Eq, r < p < q < n} = so{n — r). 

Therefore, according to the observations of Section 2, the multidimensional 
Veselova system can be treated as a generalized Chaplygin system on the prin- 
cipal bundle 

SO{n -r) — ^ SO{n) 

I TT , (5.15) 

V{r,n) = SO{n ~ r)\SO{n) 

where V{r, n) is the Stiefel variety. It can be regarded as the variety of ordered 
sets of r orthogonal unit vectors ei, . . . , in R" (C"), or, equivalently, the set 
of r X n matrices X = (ci ■ ■ ■ Cr) satisfying X^X = 1^, where 1^ is the r x r unit 
matrix. Thus V{r,n) is a smooth variety of dimension N = rn — r{r + l)/2 
(see e.g., (23; ) . and the components of X are redundant coordinates on it. 

The nonholonomic right-invariant distribution Dr is orthogonal to the leaf 
of the action of SO{n — r) with respect to the bi-invariant metric on SO{n). 

The tangent bundle TV{r, n) is the set of pairs X, X subject to the con- 
straints 

X^X^lr, X^X + X^X ^0, (5.16) 
which give r{r + 1) independent scalar constraints. 

Lemma 5.3 The momentum map $ : TV{r,n) —> so{n) is given by 
Lo = X) = xx'^ - xx'^ + ^X[X'^X - X'^X]X^ 

= eiAeiH h A -I- ^ ^ [(e^, e/s) - (ea, e^)] A e^. (5.17) 

l<Q</3<r 

Indeed, $(<%", X)'^ = -$(<%", X) and therefore ^{X, X) e so{n). Taking into 
account constraints H5.16|l . we obtain 

X)X ^X - XX^X - \^X{X^X - X'^X) , 
which implies the Poisson equations for e^, 

X^-ujX. (5.18) 
On the other hand, putting X — ~uoX into ^{X,X), we get 

$ = ljT + Tlj — TljT = prp (cj), F = ei ei -|- • • • + e^. 
Hence A") G Vr and formula (|5.17l) describes the momentum mapping. 
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The reduced Lagrangian L{X,X) takes the form 

L ^ ^{I^X,X),<1>{X,X)) = -itr (l<^>{X,X) o ^{X,X)^ . 
Then we introduce the r x n momentum matrix 

r.,s=dLiX,X)/dX,s. (5.19) 

Since the Lagrangian is degenerate in the redundant velocities Xig, from this 
relation one cannot express X in terms oi X,V uniquely. On the other hand, the 
cotangent bundle T*V{r,n) can be realized as the set of pairs XjV satisfying 
the constraints 

X'^X^Ir, X^V + V^X ^0. (5.20) 

The corresponding symplectic structure ]7 on T*V{r,n) is just the restriction 
of the canonical 2-form on the ambient space M^"'' = {XjV), 

n r 
i=l s=l 

Under conditions H5.2()(l . relation (|5.19(l can be uniquely inverted, i.e., one gets 
X — X{X,V) (for r = 1 see the section below). 

Next, according to the definition of the reduced momentum map in H2.11() 
and in view of (|5.17|l . (|5.19l) . the map $* : T^V{r,n) so*{n) has the form 

<i>*{X, V) = Iu\-D, = XV^ - VX'^ + ]^X[X'^V - V^X]X^. (5.21) 

It establishes a bijection between the linear subspace I?,- C so*{n) — {M} and 
the cotangent space T^V{r, n). 

Theorem 5.4 The reduced LR system on T*V{r, n) is the restriction of the 
following system on the space {X,V), 

X = -~Lu{X, r)X, V = -uj{X, V) V, (5.22) 

where uj{X,V) = ^X,X{X,V)). 

Proof. Substituting the expression H5.21|l into the reduced momentum equation 
H2.12|l . differentiating its left hand side, then taking into account the Poisson 
equations H5.18|l and the conditions (|5.20() . we obtain 

XV^ _ 'PX'^ _ Xi'P'^X - V'^LoX)X^ = XV^uj + iuVX^. 

Multiplying both sides of this relation from the left by X^ and from the right 
by X, then using again the conditions H5.20() . we arrive at the equation 

-X'^V ^ X^ujV, 

which implies that V = —ujP, i.e., the second equation in (|5.22|l . The first 
equation in this system is just a repetition of H5.18() . The theorem is proved. 

According to H5.22|l . apart from the energy integral, the reduced flow on 
T*V{r,n) possesses matrix momentum integral V"^V. 

Notice that the form of equations (|5.22l) is similar to those describing geodesic 
flows on Stiefel and Grassmann varieties (see [71^]) and on other homogeneous 
spaces f |46l 1141 111] '). However, our system is not Hamiltonian with respect to 
the symplectic structure Q. 



24 



Reduced invariant measure. The phase space {M,ei, . . . , e^) of the LR 
system (|5.6ll . H5.7II . (|5.8|l has the structure of the dual to the semi-direct Lie 
algebra product 

so{n) X ( R" X ■ ■ ■ X ]R"j 
r times 

and carries the corresponding Lie-Poisson structure {•, -js- (For r = 1 this is 
just the Lie-Poisson bracket on the coalgebra e*{n).) This Poisson structure is 
degenerate, and the subvariety C- C {M,ei, . . . , e^) defined by the constraints 
H5.9|) and the conditions 

4>ki ^ iek,ei) = Ski, fc, / = r 

is its symplectic leaf of dimension 2N — 2 dim V{r,n) = rn — r{r + l)/2: the 
restriction of {•, - js onto Or is nondegenerate. The variety Or is a (generally 
singular) orbit of coadjoint action of the semi-direct group product SO{n) x 
(M" X • • • X R"). 

Let S be the corresponding symplectic structure on Or- By construction, 
the extended momentum map 

: T*V{r,n) ^ Or, {X,V) ^ {ei, . . . ,er,<S>*{X,V)) 

preserves the Poisson structure, hence it is a symplectomorphism: the 2-form 
E passes to the symplectic structure J7 on T*V{r,n). Thus, , as a volume 
form on Or, transforms to the canonical volume form on T*V{r,n). 

As we know from Theorem 3.3, a reduced LR system always has an invariant 
measure. Using the above property of $*, the measure in our example can be 
written explicitly. 

Theorem 5.5 The reduced LR flow on T*V{r,n) has invariant measure 

l/Vdet(X|pJI2^, N = dim V{r,n) = rn-r{r + l)/2. 

Notice that the density of this measure coincides with that of the invariant 
measure Q of the LR system (|5.6f) - (|5.8f) in the coordinates {M, ei, . . . , e^). In 
particular, for the special inertia tensor H5.11|l the density is the same as in 

Sketch of proof of Theorem \b.5\ Let 

i/'i(A^,e), . . . ,'(/'m(Al,e), m= {n - r){n - r - l)/2 

be any independent linear combinations of the constraint functions defined by 
1)5. 9() . Then, at points of Or C {M,ei, . . . , e^), 

dMuA • • • A dMn-i,n A den A • • • A de™ 

= C(ei,...,e^)-dV'iA---AdV'mAZ'^ Y[ hd(j)ki, (5.23) 

l<k<l<r 

where, as above, (l)ki — (e^, e;) and ^(ei, . . . , e^) is a certain nonvanishing func- 
tion. Let X^.^ , Xe-^ be the basis vectors in the phase space of the LR system. 
The function ^ can be found by inserting the polyvector 

Xmi2 a • • • a -^A1„_i,„ a Xeii A • • • A Xe„r 
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into the left and right hand sides of H5.23|l and taking into account the relations 
^i^M^J,XMpg) = {Mtj,Mpq}, S{XM,j, Xej,^) = {Mtj,epk}, 

^(^epfc,^e,J {epk,eql} = 0. 

One can always choose such a basis of functions tpk (Al , e) that becomes 
a constant on the whole orbit Or- Moreover, for this basis, the time derivative 
with respect to the flow H5.6|l - H5.8|l has the form ipk = J2T=i ^ksi^s where the 
functions >Cks satisfy the conditions >Ckk = on Or- 

Now let £y = diy+i^^ d denote the Lie derivative in the space {M,ei, . . . , e^) 
with respect to this flow, iy being the interior product corresponding to the flow. 
Since the functions (jj^i are its generic first integrals, one has 

Cydc^ki=d{<i)ki) = 0- (5.24) 

On the other hand, the functions V'fe £^re particular integrals of the flow. Then 
for the chosen basis of such functions and for any k, 

m 

Cydipk = d{i>k) = ^(^ffcsdV's + ipsdxks) with Xkk = on Or - 

s=l 

Hence, at points of Or we have 

dV'i A---A(£i,dV'fc)A---AdVm H Ad0fe,=O. (5.25) 

l<fe</<r 

Now, combining relations (|5.10() and H5.23|l with =const, we get 

CyQ = Cy ( 1/Vdet(l|cj d^Ai A • • • A d^Am A Y[ Ad^bkij^O. 

This, together with f05|l . implies that £y{l/ y/det{l\vj U^) = 0; 

therefore the volume form under the Lie derivative is an integral invariant. In 
view of the symplectomorphism between Or and T*V{r, n), we can replace the 
volume form by which yields the statement of the theorem. 

Reducibility of the system (|5.22|l to the Hamiltonian form via a time rescal- 
ing for an arbitrary rank r > 1 and an arbitrary inertia tensor 2 is an open 
problem. 



6 Rank 1 case and integrable geodesic flow on 

gn-l 

Now we concentrate on the case r = 1 given by the original condition H5.2I) and 
again assume that the inertia tensor has the form (|5.11() . The variety V{l,n) 
can be realized as the unit sphere 5'""^ in M" = (gi, . . . , (/„), 

S-'={ql + --- + q^ = l}, 

where we set q = ei, and the momentum map H5.17|l is reduced to 

LU = <P{q,q) = qAq- (6.1) 
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Therefore, for the solution (ei(t), uj{t) — ei{t) Aei(i)) of the system (|5.6ll (|5.8|l . 
q{t) = ei{t) is a motion of the reduced system on the sphere S"""^. 

For the analysis of the reduced system we can use Theorems 15.41 15.51 of 
the previous section. However, for our purposes we shall use the reduction 
procedure described in Proposition 2.1. 

Under the condition H5.11() on the inertia tensor and in view of H6.1() . the 
reduced Lagrangian L{q, q) and the right hand side of the Lagrange-d'Alembert 
equation (|2.1U|) take the form 

L = j^[{Aq,q){Aq,q)~{Aq,qr], (6.2) 



1 

detA 

1 1 



(X$((7, q), prg-i^g[$(9, q), $(q, 0]) = -rr^iM A Aq, prg-i[,g C A g) 



^^^^{q,Aq)iAq,0 - ^^{q, Aq){Aq, = ^{q,q,0- (6.3) 
Here we used the relation prg-i|,g ^ A g = ^ A q for any admissible vector ^ = 

(ei,...,C„)^eT,5"-i. 

As in Section 4, the reduction of the LR system H5.4fl onto T*S"~^ can 
explicitly be written in terms of local coordinates qi, . . . , q-n-i on 5""^ and the 
corresponding momenta. 

As an alternative, below we shall keep using the redundant coordinates qi 
and velocities gi, in which the Lagrange equations have the form 

dL d dL . 

"a 17"a^ = '^» + ^9*' i = l,...,n, (6.4) 

oqi at oqi 

9* 1 1 

where A is a Lagrange multiplier. We now want to represent the reduced 
LR system on T*S'^~^ as a restriction of a system on the Euclidean space 
K.^" = Note that L{q,q) is degenerate in the redundant velocities q, 

hence they cannot be expressed uniquely in terms of the redundant momenta 

= II " diL^'^'^'^)^''^^ - di:4(^'^'^^^''^' ■ ^'-'^ 

In this case one can apply the Dirac formalism for Hamiltonian systems with 
constraints in the phase space (see, e.g., j2]|21Eni). Namely, from H6.5|l we find 
that {q,p) = 0, hence the cotangent bundle T*S^~^ is realized as a subvariety 
of K^" = (q,p) defined by the constraints 

01 = {q, q) = 1, 02 = {q,p) = 0. 
Under these conditions, relations IjCSfl can be uniquely inverted to yield 

det^ 



(g, Aq) 



[A-'p~ip,A-'q)q] . (6.6) 



Next, we note that dL/dqi ~ tt^. Then, from H6.4|l we obtain p = — Ag, and 
from the condition {q^p) + [q^p) ~ 0, 

p = - Ag, A = det ^ fa-^'M - _ 

{q,Aq) 
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The system (|6.6|l . (|6.7|l on r*5'" ^ coincides with the restriction of the 
following system on M^" ~ {q,p}' 

Tr.t{q,p) ^TTi{q,q{q,p)), H = ^detA ^^f^^j^J^^ ■ 

Here {F, G}* denotes the Dirac bracket on R^", 

{F,cj,^}{G,(t,2}-{F,cj,2}{G,M 



{F, G}, = {F, G} 



{01, 02} 



and {•, •} is the standard Poisson bracket on R^". The latter system has the 
explicit vector form 



q = 



det A 



A-^p- 



{p,A-'q) ^ 

iq,q) 



(p,^ ^p){q,q)- {p,q){q,A ^p) 
p = - det A — -— q . 

iq,M)[q,qr 

The bracket {•, is degenerate and 0i, 02 are its Casimir functions. 
Notice that from H6.1|l and (|6.6|l we get 



(6.8) 



uj ^ q /\ 



dot A 

(q, M) 



, jp.A-^q) ^ 

^ p 7 ^ q 



Then equations (|6.8|l can also be obtained directly from Theorem 5.4 by setting 
r=l,X = q,V=p. 

Finally, from Theorems 5.2 and 5.5 we get the following corollary. 

Corollary 6.1 The reduced LR system il6'.6)) . j6'.7| ) on T*S"^^ possesses an 
invariant measure 

(A<z,g)-("-2)/2^^ a = 

where a is the volume 2(n — l)-form and fl is the restriction of the canonical 
symplectic form dpi A dqi + • • • + dpn A dqn onto T*S"^^. 

In particular, for the reduction of the Veselova LR system (|4.2|l onto T*S^, 
the density of its invariant measure is proportional to {q, Aq), as was al- 
ready claimed in Section 4. 



Reducibility. As follows from Corollary 16.11 item 1) of Theorem 13.51 and 
the fact that the dimension of the reduced configuration manifold equals n— 1, if 
our reduced LR system on T*S'"~^ were transformable to a Hamiltonian form 
by a time reparameterization, then the corresponding reducing multiplier J\f 
should be proportional to 1/ a/ {q, Aq). 

Although Chaplygin's reducibility theorem does not admit a straightfor- 
ward multidimensional generalization, i.e., item 1) of Theorem 13.51 cannot be 
inverted, remarkably, for our reduced LR system on T*S"'~^ the inverse state- 
ment becomes applicable. 
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Theorem 6.2 1). Under the time substitution dr — ^det A/{Aq, q) dt and 
an appropriate change of momenta, the reduced LR system \6.4\j or ^6'.6]) . 
{6. 7| ), becomes a Hamiltonian system describing a geodesic flow on S"~^ 
with the following Lagrangian obtained from |6'.ij|] 



L*{q,dq/dT) = ^{q,Aqr' 



(6.9) 



2). The latter system is algebraic completely integrable for any dimension n. 
In the spheroconic coordinates Ai, . . . , A„_i on 5"^^ such that 

^= ^'''!f--:^'\':"-'\ h-Ar\ (6.10) 

the Lagrangian L*(q,dq/dT) takes the Stdckel form and the evolution of 
Afc is described by the Abel-Jacobi quadratures 



2/rIAO 2v/i?(A„-i) 
k — 1, - ■ ■ , n — 1, 

where 

i?(A) = -(A - /i) • • • (A - /„)A(A - C2) • • • (A - c„,i), (6.12) 

h = L* being the energy constant and C2, ■ ■ ■ , c„_i being other constants 
of motion (we set Ci = 0). For generic values of these constants the 
corresponding invariant manifolds are {n — 1)- dimensional tori. 



We start with the proof of item 2) of Theorem 16. 21 which is quite standard. 
Using the Jacobi identities, 

^ '0 forO<rn<7V-l, 



for any distinct pi , . . . , pAr, = = <^ 



i^iUi^siPi - Ps) U forrn = iV-l, 
in the spheroconic coordinates we have 



1 Ai---A„_i Yls^ki^k-Xs) (d 



4 /i • • • /„ ^ 

^ " fc=l 



(A/c - /i) • • • (Afe - /„) Afc Vc^'T 



(Ag,g)^(/-^g,g)^ \---Y' - (6-14) 
Then the reduced Lagrangian L*{q,dq/dT) in H6.9fl takes form 

1'^^ n.^fc(A,-A.) /d 

8fr[ (Afe-/i)---(Afe-/„)Afc \dT 

As a result, the corresponding Hamihonian written in terms of 

dL* 



Afe, Pfe — 



5(dAfe/dT) 
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is of Stackel type (see e.g., jl]), which leads to the quadratures H6.11() and proves 
the integrabihty of the system. 

The proof of item 1) of Theorem 16.21 is based on a relation between the 
reduced LR system and the Neumann system on S"^^ and will be given at the 
end of this section. 



Reduction to the Neumann system. It appears that Theorem 14.11 
relating the Veselova LR system and the classical Neumann system has the 
following multidimensional generalization. Introduce another new time ti by 
formula 

dr, = r'dt, = ^^^^^^^^^Ilpldt, (6.15) 

and let ' denotes the derivation in the new time. 

Theorem 6.3 Under the time substitution f6.15\) . the solutions q(t) of the 
reduced multidimensional Veselova system on S""^^ transform to solutions of 
the integrable Neumann problem with potential U{q) — ^{A''^q,q), 

q"^-jq + Xq, = ^ > (6-16) 

corresponding to the zero value of the integral 

Fo{q,q') = {Aq',q'){Aq,q) - {Aq,q'f - {Aq^q) (6.17) 

and vise versa. 

The proof we shall present here is similar to that of "three-dimensional" 
Theorem 4.1, which was given in |48j . 

Proof of Theorem 16 .HI Let w = g A g as above, and set P = q h q' . Then the 
energy integral of the reduced Veselova system and the integral (|6.17|) of the 
Neumann system can be written as 

E{q,q) = ^{Ju;,Lo), F{q,q') = detA{IP,P) - {Aq,q). 

The time reparameterization (|6.15|l induces a bijection between invariant 
submanifolds ^ {E ^ h} d TS"-^{q, q} and jPg = {F = 0} C TS'^-^iq, q'}. 
Indeed, on £h we have 



_i l2hdctA 



dt = iihdTi, =Y"(I7^- (^-^^^ 

Then the point (<?,(?) € £h corresponds to {q,q'), q' = /i/ig, and the equation 
{Iuj,uj)/2 ~ h corresponds to the relation 



2^2^ ' 2 iAq,q) ^ ' 



Therefore F = det A{IP, P) - (Ag, q) = 0, and {q, q') G J^o- 

Next, note that equations (|5.7() with uj = q A q are equivalent to 

{Iclj-q)Aq+{Iuj-q)Aq = 0. (6.19) 
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In view of time reparameterization (|6.18() we have that P = iih^ and 

dP dP d 2dLLl ^ d 2\ Ir. nr^\ 

rf7 ^ dF'^'' ^ = A^.^ + -^J^^^n)^ ■ (6-20) 

Now we apply the inertia operator H5.11|l to both sides of this relation, then 
multiply the result by the vector q, and finally take the wedge product with q. 
As a result, due to (|6.18|) . we get 

Using H6.19|l . we transform H6.21|l into 

2ft. det ^ {IP' ■ q) A q = -{Aq,q){Iuj ■ q) A q + {Aq,q){Iuj ■ q) A q. (6.22) 
The right hand side of H6.22|l is of the form E A q, where 
S = {Aq, q)XLJ ■ q — {Aq, q)Xuj ■ q 

= -^^{Aq,q){Aq®Aq-Aq®Aq)-q 

-—^{Aq, q){Aq ® Aq - Aq ® Aq) ■ q ^ -2hAq. (6.23) 
det Ji 

For the last equality in l|6.23|l we used the identity 

2h= {I{qAq),qAq) = -^^{Aq A Aq, q A q) = -^^{Aq, q){Aq, q) ~ {Aq, qf . 
Hence, and yield 

{IP' ■q)Aq^ ^^qAAq, P = q A q' . (6.24) 

In view of the constraint {q, q) — 1, this is equivalent to equations (|6.16|) . 

Thus we proved that if q{t) is a solution of the reduced multidimensional 
Veselova system laying on Eh, i.e., q{t) satisfies (|6.19|l . then q{Ti) is a solution 
of the Neumann system (I6.16|) laying on J-'q. 

Conversely, starting from H6.24|l and repeating calculations in the inverse 
direction, we arrive at (|6.19|) . The theorem is proved. 

It is known (see e.g., |35l I40L that the Neumann system on S"~^ pos- 
sesses the following family of quadratic first integrals: 

p2 "2 

^Si^^ («^^^) 

l<i<j<ri ^ ' 1=1 

and that the evolution of the spheroconic coordinates defined by (|6.10() is 
described by equations 

where TZ{X) is the following polynomial of degree 2n — 1: 

n = -$2(A)^(A), $(A) = (A - /i) • • • (A - /„) . 
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Next, as follows from l|6.25|l . for the trajectories (?(ti) corresponding to the 
zero value of the integral H6.17(l . we have ^{0) — 0. Hence in this case, the 
polynomial TZ{X) has the same form as H6.12|) . that is 

7^(A) - -(A - /i) • • • (A - /„)A(A - C2) • • • (A - c„_i) . (6.27) 

Now, comparing equations H6.26|l with the quadratures H6.11|l . we arrive at 
the following proposition. 

Proposition 6.4 Under the time rescaling dri — \/2hdT the solution <?(ti) of 
the Neumann system \(i.lf)]) lying on Tq — {-Fq = 0} transforms to a solution 
q{T) of the geodesic flow on S""^^ described by the Lagrangian L* in 1^6. y\) and 
having the energy constant h, and vise versa. 

Now combining Theorem lG . 31 and Proposition l6.4l we finally obtain the proof 
of item 1) of Theorem 16.21 



7 Reconstructed motion on the distribution D 

Here we study the integrability of the original (unreduced) LR system on the 
right-invariant distribution D C TSO{n) of dimension (n — 1) + n{n — l)/2, 
which is specified by constraints H5.2|l and the left-invariant metric defined by 
llOTTl . 

In the Hamiltonian case, the integrability of the reduced system generally 
implies a non-commutative integrability of the original system, namely its phase 
space is foliated by invariant isotropic tori with quasi-periodic dynamics. In 
our nonholonomic case one has to solve the reconstruction problem: find all 
trajectories {g{t),g{t)) in D that under the 5'0(n— l)-reduction tt : D TS^^^^ 
are projected to the given trajectory {q{t),q{t)) in TS^^~^. (In particular, for the 
Fedorov-Kozlov integrable case of the multidimensional nonholonomic Suslov 
problem, the reconstruction problem was studied in |50[ 

Since SO{n — 1) is a symmetry group of the LR system on D, and the 
reduced motion on T5"^^ occurs on (n — 1) -dimensional generic invariant tori, 
it is natural to expect that the reconstructed motion {g{t), g(t)) is quasi-periodic 
over {p-\-n— l)-dimensional tori, where p does not exceed the dimension of the 
maximal commutative subgroup of SO{n — 1), that is p < rank SO{n — 1) = 
[^] (see m)- 

As we shall see below, for our case this is not true. In fact, the relation 
between the reduced LR system and the Neumann system described by Theo- 
rem 16.31 enables us to reconstruct the motion on D exactly. For this purpose 
we also shall make use of the remarkable correspondence between the Neu- 
mann system and the geodesic fiow on a quadric. Namely, consider a family of 
(n — 1) -dimensional confocal quadrics in R" — {Xi, . . . , X„), 

[a - Ai a - An ) 

where Ai, . . . , An are distinct numbers. 

Theorem 7.1 f|35p. Let X{s) be a geodesic on the quadric Q{0), s being a 
natural parameter. Then under the time rescaling 



{dX/ds,A~^dX/ds) ^ 
^'=V iX,A-^X) ^^-2) 
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the unit normal vector q{Ti) = A^^X/\A^^X\ at the point X E Q is a solution 
to the Neumann system 16.16]) corresponding to the zero value of the integral 
Fo{q,q') in \6.16]) and vise versa. 

Next, it is well known that the problem of geodesies on a quadric Q(0) 
is completely integrable and that the qualitative behavior of the geodesies is 
described by the Chasles theorem (see e.g., j^SlQUI)- the tangent line 



of a geodesic X{s) on (5(0) is also tangent to a fixed set of confocal quadrics 
Q{a2), ■ . ■ , Q(ckn~i) C M", where a2, . . . , a„-i are parameters playing the role 
of constants of motion (we set ai — 0). Now let be the normal vector 
of the quadric Q{(y^k) at the contact point pk = (■ C] Q{ak)- Then another 
classical theorem says that the normal vectors ni, . . . ,n„_i, together with the 
unit tangent vector 7 = dX/ds, form an orthogonal basis in R". 

The integrability of the geodesic flow is also related to the following prop- 
erties found by Moser in |4()j . 

Proposition 7.2 1). Let x be the position vector of any point on the line is, 
which is tangent to the geodesic X(s). Then in the new parametrization 
si such that ds — ~{X, A^'^X) dsi the evolution of the line is described 
by Lax-type equation in n x n matrix form 



where Tl^ = Id — (7, 7) "'^7 7 is the projection onto the orthogonal 
complement of 7 in R" . 

2). The conserved eigenvalues of C are given by the parameters ai = 0, a2, . ■ . , ctn-i 
of the confocal quadrics and by an extra zero. The corresponding eigen- 
vectors of C are parallel to the normal vectors ni = q, . . . , n„_i, and to 
1- 

Now we are ready to describe generic solutions of the original LR system 
on C TSO{n). Let q{Ti) be the solution of the Neumann system Ht).15|l with 
-fo(<7;9') = 0, which is associated to a solution {q{t),p{t)) of the reduced LR 
system as described by Theorem 16.31 Let 



X = iq,Aq)-'/^Aq{s), m - g(s), . . . , n„_i(s), lis)^— (7.5) 



be the corresponding geodesic on Q{0) in the parametrization s given by (|7.2fl 
and the unit eigenvectors of £ associated to this geodesic. Note that according 
to l|().15|l and (|7.2I) . one can treat s as a known function of the original time t. 
Then the following reconstruction theorem holds. 

Theorem 7.3 For a solution {q{t),q{t)) of the reduced LR system on T S"^^, 
a solution g(t) of the original system on the distribution D <zT SO{n) is given 
by the orthogonal frame formed by the unit vectors 



Is = {X{s) + CT dXjds I (7 e M} 




B = A~^x® A-^ 



£ = (A - a; (g) a;)n. 
^7- A"^7(g) A^^a;, 



(7.3) 



(7.4) 



ei = q(t), 62 = n2(s(t)), . . . , 



e„-i = n„_i(s(i)), 



e„ = l[.s{t)). 
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The other solutions {g{t), g{t)) that are projected onto the same trajectory {q{t), q{t)) 
have the same ei, while the rest of the frame is obtained by the orthogonal trans- 
formations 

(ezW-'-enW) = (n2W---n„_i(i)7W)^, (7-6) 
where d\ is a constant matrix ranging over the group SO{n — 1). 

From Theorems 17.31 and the integrabiUty properties of the Neumann system 
on r*S'"-i (Theorem ESI), we conclude that the phase space D C T SO{n) of 
the muhidimensional Veselova LR system with the left-invariant metric defined 
by H5.11|l is almost everywhere foliated by {n — l)-dimensional invariant tori, 
on which the motion is straight-line but not uniform. 

This also implies that, apart from the pull-back of the n — 1 integrals of the 
Neumann system, the LR system possesses [n — l)(n — 2)/2 generally indepen- 
dent integrals on D. In particular, as follows from the nonholonomic momentum 
equations H2.9|l with Wfc = ei A Cfc and the special form of the inertia tensor, 
the LR system on D has linear integrals 

Ik = (Al,ei A efc) ee -i— -(AwA,ei A e^), k^2,...,n, 
det A 

of which n — 2 ones are independent, since ^2 + ' ' ' + = P^- 

Proof of Theorem \7 .'6\ As follows from ProDOsition l7.2l for the geodesic motion 
on Q(0), the unit normal vectors rii, . . . , n^-i and 7 satisfy the equations 

-^nfc = -Snfc, --^7 = -B7, /c = 1, . . . , n - 1, 
asi asi 

ds = -vdsi, ly = (X,A^'^X) , 

which can be regarded as kinematic (Poisson) equations with the "angular 
velocity" matrix B. 

Next, from Theorem 17. II we have 

X = V^Aq and ^^^=^A^ + ^Aq. (7.7) 

ds ds ds 

Now let us choose x = X{s) in the expression 1)7. 4|l for B. Substituting (|7.7(l 
into this expression we find that B — v q A dq/ds. Then, in the original param- 
eterization s, the above Poisson equations take the simple form 

dnk/ds — ~{q A dq/ds)nk, d^y /ds = — {q A dq/ds)j . 

Changing here the time parameter s to t and taking into account relation H6.1|l . 
we finally obtain 

hk = -uj{t)nk, fc = l,...,n-l, 7 = -tj(t)7, (7.8) 

where w G I? C so{ji) is the admissible angular velocity of the n-dimensional 
body. This implies that the orthogonal frame {ni{t), . . . , n„_i(t), 7(t)) gives a 
solution of the LR system on D. 

Note that the vectors of the frames that are obtained by the orthogonal 
transformations (|7.6() also satisfy the Poisson equations H7.8|l and therefore also 
give solutions of the LR system on D. Since the fiber of the map tt : D ^ TS'^~^ 
is the group SO{n — 1), there are no other solutions on D that are projected 
onto the same trajectory {q{t),q{t)). The theorem is proved. 
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The explicit solution for the frame {q{t),nk{t),^{t)). In order to 
find explicit expressions for the components of Wk and 7, following Jacobi 
we first introduce ellipsoidal coordinates j^i, . . . , Vn-i on Q{0) according to the 
formulas 

Matching these with the expressions H6.10|l for qi in terms of the spheroconic 
coordinates Ai, . . . , A„_i on 5*"^^ and taking into account (|7.5|l and (|6.14|l . we 
find that, up to a permutation of indices, 

Vk = AjT^ fc = 1, . . . , n - 1. 

Using this property one can also prove that the nonzero parameters 122, • ■ • > Q;„_i 
of the confocal quadrics in the Chasles theorem are just the inverses of the con- 
stants of motion C2, . . . , c„_i in the invariant polynomial H6.27|l . 

7^(A) - -(A - /i) • • • (A - /„)A(A - C2) • • • (A - c„_i) . 

As a result, making use of the definition of the vectors x\.k^ 7, one can express 
their components in terms of pairs (Ai, a/7?,(Ai)), . . . , (A„_i, •\/7?.(A„_i)) and 
C2, . . . , c„_i. After some calculations we get 



V£/(A,/,) vi7(a:o) 



where 



j = l,...,n, fc = 2,...,n— 1, 
^'(r) = (r - /i) • • • (r - /„) , ip{r) = r(r - C2) • • • (r - c„_i) , 



C/(A,r) = (r-Ai)---(r-A„_0, = . 

n(A.-A,) 

The evolution of A-coordinates in the time r is described by the quadratures 

The squares of expressions (|7.9f) are symmetric algebraic functions of n — 1 
points (wkjXk) on the hyperelliptic curve C — {w^ = 7?.(A)} of genus n — 1. 
Then, by using the classical algebraic geometrical methods (see, e.g., 01201); 
the components of g, n2, . . . , n„_i, 7, and the function -^Ai • • • A„_i can be rep- 
resented as quotients of theta-functions with half-integer theta-characteristics 
associated to C, whose arguments depend linearly on ti. 

Finally, the dependence of t in ti is obtained by the integration of H6.15|l . 
which, in view of H().14|l . leads to the simple quadrature 

^ V '^i('''i) • • • ^n-iiri) dri + const. 
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Conclusion 



In this paper we considered LR systems on compact Lie groups and showed that 
their reductions to homogeneous spaces always possess an invariant measure. 
We calculated the density of this measure explicitly in case of the Stiefel varieties 
V{r,n) = SO{n)/SO{n — r). It turned out that for r = 1 and the special 
inertia tensor on so{n) , the reduced flow is transformed to an integrable geodesic 
flow on S"'~^ via the time rescaling prescribed by the density of the invariant 
measure and Theorem lr>.5l Moreover, in this case the unreduced flow on the 
right-invariant distribution D G T SO{n) is also integrable. 

Such a behavior of a multidimensional nonholonomic system is exceptional, 
which may follow from the rich underlying geometry coming from the Chasles 
theorem and the Jacobi problem on geodesies on an ellipsoid. (The latter is 
closely related to integrability of various problems in mechanics. 

In this connection the following questions arise: are there other inertia ten- 
sors of LR systems on SO{n) for which the above properties hold and how wide 
is the class of such tensors? Can reduced flows on V{r, n), r > 1 be transformed 
to the Hamiltonian form (with respect to the canonical symplectic structure J7) 
in the same manner? Are the reduction to a Hamiltonian form and integrability 
still possible in case of nonhomogeneous right-invariant constraints on SO{n) 
(similar to what takes place for the classical case n = 3)? 

A part of our analysis can be extended to LR systems on noncompact Lie 
groups and their reductions. It would be interesting to study meaningful ex- 
amples of such systems. 
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